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ABSTRACT 

The whole of the Taurus region (a total area of 52 sq. deg.) has been observed by the Her- 

schel SPIRE and PACS instruments at wavelengths of 70, 160, 250, 350 and 500 /im as part 

of the Herschel Gould Belt Survey. In this paper we present the first results from the part 

of the Taurus region that includes the Barnard 18 and L1536 clouds. A new source-finding 

routine, the Cardiff Source-finding AlgoRithm (CSAR), is introduced, which is loosely based 

on CLUMPFIND, but that also generates a structure tree, or dendrogram, which can be used to 

isj> ' interpret hierarchical clump structure in a complex region. Sources were extracted from the 

'V^ I data using the hierarchical version of CSAR and plotted on a mass-size diagram. We found a 

d , hierarchy of objects with sizes in the range 0.024-2.7 pc. Previous studies showed that gravi- 

tationally bound prestellar cores and unbound starless clumps appeared in different places on 
the mass-size diagram. However, it was unclear whether this was due to a lack of instrumental 
dynamic range or whether they were actually two distinct populations. The excellent sensi- 
tivity of Herschel shows that our sources fill the gap in the mass-size plane between starless 
and pre-stellar cores, and gives the first clear supporting observational evidence for the the- 
ory that unbound clumps and (gravitationally bound) prestellar cores are all part of the same 
population, and hence presumably part of the same evolutionary sequence (c.f. Simpson et al. 
2011). 

Key words: stars: formation — ISM: dust — infrared: ISM — ISM: individual (Taurus, 
TMC2, Barnard 18) 



1 INTRODUCTION 

* //eric/je/ is an ES A space observatory with science instruments provided r,,, ^,,tii.- i / ii j- i\i r i i 

. . . . , . , . . The Gould Belt is a large (~ 1 kpc diameter) loop oi molecular 

by European-led Pnncipal Investigator consortia and with important paitic- , „„ ■ ■ ■ ■ ,■ nr,o ^ 

. ^. „ xTACA clouds and OB associations that IS inclined at ~ 20 to the Galac- 

ipation trom NASA. i ii -— i i — ni —t 

t Email: jmkirk@uclan .ac.uk tic Plane jHerschelll 18471 ; iGoulal 1 8791) . The Gould Belt is impor- 
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Figur e 1. A coverage map of the Herschel Gould Belt Survey jAndre et al.1 
1201011 of the Taurus star fonnation region. The greysca le and contours 
show visual extinction calculated from the 2MASS survey iSchneider et alj 
I2OI 11) . Plotted over this are the foot-prints and names of the Herschel fields. 
The regions presented in this paper are outlined by the thick lined boxes. 
Note that this image covers ~ 10° X 10°. A single A^ = 3 contour is 
plotted. 



Table 1. Multi-wavelength properties of the Herschel maps. The first line 
lists the native instrumental FWHM. The second fine list the median ab- 
solute monochromatic intensities towards the S1-S2 fields as derived from 
IRAS/Planck observations. The third line lists the Icr rms levels as mea- 
sured from the maps after they have been convolved/resampled to the 
500/xm resolution and pixel-grid. 



Property 


500/im 


350^m 


Wavelength 
250Afm 


160^m 


70^m 


Beam FWHM (") 
Median (MJy/sr) 
lo- (MJy/sr) 


36 
10.7 
0.8 


25 
23.4 
1.5 


18 
38.8 

2.2 


12 X 16 
60.0 
2.9 


6 X 12 
3.97 
3.8 



tinuum and clumps seen by molecular line tracers such as CO 
lie in very different regions of this plane. There are two possi- 
ble explanations for this difference: either selection effects are at 
play, or there are two fundamentally different populations of ob- 
ject, i.e., one of low densit y clumps that rarely form stars (c.f. 
IWard-Thompson et aLluOlOf) . and one of higher density cores that 
are the pr ecursors to star formation, no rmally referred to as prestel- 
lar cores (IWard-Thompson et al.ll994ii2007) . In this paper we will 
use the first look Taurus dataset to test Herschel's ability to populate 
this mass-size plane and examine whether it is possible to demon- 
strate which of the two posssible explanations is more likely. 



tant for star formation studies as it contains most of the nearby 
(d < 500 pc) low-mass star formation regions, where our obser- 
vations have their highest spatial resolution and mass sensitivity. 
Thus the Gould Belt is an excellent laboratory in which to study all 
stages of low-mass star formation. For these reasons, s eparate con- 
sortia from the Herschel ( Andre et a l. 2010J), Spitzer JEvans et al.l 
1200^ Allen et al. in prep), and JCMT JWardThompson et akbOOTh 
communities have mapped, or are mapping, significant fractions of 
the Gould Belt at wavelengths from 3.6 |J.m to 850 i^m. 

This paper p resents the first Herschel Gould Belt Survey 
JAndre et alj|2010n observations of the Taurus Molecular Cloud. 
This cloud contains '--^ 2 — 4 x 10* M0 of molecular mate- 
rial d i stributed in a series of filaments and cores (IGoldsmith et al.l 
120081 : iKenvon et al. 2008 ). The observed star formation efficiency 
in Taurus is ~1% (Goldsmith et al. 2008). An analysis of Planck 
data towards Taurus found line-of-sight averaged dust tempera- 
tures of ~17K and ^14 K in the diffu s e and dense parts re spec- 
tively dPlanck Collaboration et ai]|201ll) . iRebuU et alj klQl&i used 
Spitzer observations to compile a catalogue ~350 new and pre- 
viously known YSO candidates across t he Taurus cloud. The es- 
tablished distance t o Taurus is 140 pc ( iStraizvs & Meistaslll98CI : 
iKenvon et alj 19941 : iLoTnard et al J200^ . 

Fig [T] shows an extinction map of Taurus produced from 
2MASS data (Schneider et al. 2011), overlaid with the fields in our 
Herschel survey. The Herschel fields cover the northern and south- 
em Taurus filaments (labelled N1-N3 and S1-S4 in our survey) and 
five satelhte complexes (L1489, L1517, L1544, L1551, and IRAM 
04191; not shown in Fig[T). The total area is 52 sq. deg. making 
it the largest region mapped in the Herschel Gould Belt Survey. 
In this paper, we present data for the first two fields from Taurus 
to have been observed, reduced, and analysed, namely SI and S2. 
These two fields are outlined in thick black boxes in Fig[T] 

Populations of molecular cores found by different tracers in 
star forming clouds occupy different regions of the mass-size plane. 
iMotte et al.l ( 1200 ih found that cores seen in the submillimetre con- 



2 OBSERVATIONS 
2.1 Data Reduction 



Herschel is a space telescope with a 3.5m-diameter mirror that op- 
erates in the far-infrared and submillimetre regimes (IPilbratt et alJ 
[2OIO). The observations for this paper were taken simultaneously 
with the Photodetecto r Array Camera and Spectrometer, PACS 
dPoglitsch et al.l l2010h, and the Spectral and Photometric Im ag 



13), and the Sp e 
i lGriffin et alM 



ing Receiver, SPIRE JGriffin et alj(2010l : ISwinvard e7aLll2010h us- 
ing the combined fast-scanning (60 arcsec/s) SPIRE/PACS paral- 
lel mode. These data used PACS filters centred at wavelengths 
of 70-/xm and 160-/im, with angular resolutions of 6"xl2" and 
12" X 16" respectively (accounting for the scan speed), and SPIRE 
filters centred at wavelengths of 250 /im, 350 /^m, and 500 /im, 
with angular resolutions of 18", 25" and 36" respectively. The an- 
gular resolutions are also listed on line 1 of Table[T] The final PACS 
and SPIRE angular resolutions are equivalent to spatial resolutions 
of 0.008-0.024 pc (1,700-5,000 AU) at the assumed distance to 
Taurus of 140 pc. 

The data were taken on Observing Day 275 (2010 February 
13th) and have Herschel Observation IDs of 1342190652-55. The 
total observation duration was 8 hours, 42 minutes. Flicker noise 
arising from the instrument electronics, so called 1/f noise, intro- 
duces a slowly varying component to the detector timelines. This 
was suppressed by mapping each field in two orthogonal scan direc- 
tions. The data timelines from each map were reduced separately 
and then combined at the map-making stage. 

The SPIRE data were processed with the software package 
HIPE, version 8.1 (Ott 2010) until the 'level 1' stage. The SPIRE 
observation baseline was removed using the 'destriper' module in 
HIPE and a coiTection for the relative gain of each SPIRE bolome- 
ter was applied. Maps were reconstructed using the 'naive' map- 
making algorithm. The extended source calibration was used. The 
rest of the PACS reduction and map reconstruct ion was done using 
version 16 of SCANAMORPHUS dRousselllJOlih . These maps were 
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Figure 2. Maps of the SI and S2 regions. Left: A false-colour image with red, green and blue showing the 500-/xm, 250-/im and 160-fim data respectively. 
Middle and Right: Column density (middle) and colour temperature (right) maps of the SI and S2 regions derived from the Herschel 70-/jm to 50Q-fj,ra 
and Spitzer 24-/im data iGiidel et al. 2007; see text for details). The contours show column density stalling at 5-(t with each subsequent contour level being 
1.5 times the previous level. The colour bar indicates the scale. Two large elongated filaments can be seen, along with smaller-scale structure. The northern 
filament is known as Barnard 18 or TMC 2 & 3. The southem filament is known as L1536 at its western end and reflection nebula GN 04.32 .8 at the east ern 
end. The bluish point source coincident with the reflection nebula in the left panel is the star HP Tau. The clumps within B 18 are labelled A-J JHeveill988h in 
the middle panel. The clumps in L1536 are numbered and the reflection nebula is also marked. 



found to be within 5% of the eq uivalent re construction using the 
alternative 'Madmap' mapmaker dOtnuOlOh . 



2.2 Map morphology 

The left-hand panel of Fig |2] shows the Herschel data of the S 1 
and S2 region as a false colour image. The three wavelengths used 
to make this image are 500 ^m (shown in red), 250 ^m. (shown 
in green) and 160 /im (shown in blue) convolved to the 500 pim 
resolution. The data show two approximately parallel filamentary 
structures, with the northern filament (coincident with the S2 field) 
being significantly longer than the southem filament (coincident 
with the SI field). In addition, many smaller filaments and arcs can 
be seen emanating from the two main filaments. Many warm com- 
pact sources, which appear blue, can also be seen within the data. 
The middle and right-hand panels of Fig |2] show the column den- 
sity and colour temperature maps of the same region derived from 
the data through smoothing the data to a common 36" resolution 
and fitting the resulting SED of each common pixel with a grey- 
body (/3 = 2; see Section IbTI for more details). Once again, the 
two main filaments can be seen. Objects mentioned in this paper 
are labelled on the middle panel of Fig|2] 

The northern filament has a projected length of ~ 3° (~7 pc) 
and appears roughly linear on the plane o f the sky. Thi s filament 
is known as Barnard 18 (hereafter B18: iBamard| 19191), and was 
divided into clumps labelled B18-A to B18-I bv lHev^J1988l) us- 



ing a map of ^"'CO emission; these clumps are indicated on Fig|2] 

Clump B18-E splits into two on our maps and they are labelled 

El and E2 here. B18-F is a faint feature in the Herschel data and 

is not labelled. B18-J is a new feature seen here for the first time 

and is labelled on Fig|2] The southem filament splits roughly into 

two halves. The western half is usually known as Lynds Dark Neb- 

I II — 

ula 1536, or L1536 (Lvnds 1962). The eastern half is known as the 

reflection nebula GN 04.32.8 ( Magakian 2003) . The morphology of 

this filament appears to change in Fig|2]from L1536 to GN 04.32.8. 



L1536 appears broadly similar to the clumps in the B18 filament, 
but in the region of GN 04.32.8 the filament is fainter. 

The reflection nebula GN 04.32.8 contains the star HP Tau 
which appears as a point source suiTounded by a crescent shaped 
feature in the Herschel map s, suggestive of inte raction between the 
point source and the cloud. iTorres et al.l (120091) used the VLBA to 
measure a parallax distance to HP Tau of 161.2 ± 0.9 pc. The 
interaction of the star with the cloud would imply that this part 
of the cloud is also at that distance. The plane-of-sky separation 
of the two filaments is ~2° which, at 161 pc, would equate to 
about 5.5 pc. There is some suggestion in CO surveys, however, 
that the filaments are not at t he same distance dStoiimirovic et al.l 
120071 ; iNaravanan et al.l 1200 8). Given that uncertainty we use, as 
stated earlier, the establis hed distance of 140 pc for th e entire 
cloud JStraizvs & Meistas|[r980l ; lKenvon et alJll994l ; lLoinard et al.1 
^008). If the southem filament is at ~ 160 pc and the northem fil- 
ament is closer to the canonical ^ 140 pc then the two filaments 
could be separated by up to ^--^20 pc along the line of sight. 

Comparing the left-hand panel of Fig |2] with Herschel ob- 
servations of some other regions in the Gould Belt, interesting 
differences appear. For example, compared to the Aquila region 
(Andre et al. 2010), the region displayed in Fig|2]shows relatively 
little change in colour, and thus temperature, between different 
parts of the cloud. The coldest, densest parts appear slightly red- 
der than the average, with only isolated points of blu e (i.e., warm 
YSOs). In contrast, the equivalent image of Aquila JAndre et alj 
[20 10) shows far more variation, with large areas appearing red and 
the warm W40 nebula showing up clearly in blue. This difference is 
not an artefact of the false-colour imaging; it represents a real dif- 
ference between the Aquila and Taurus star-forming regions . There 
is far more variation in temp erature across the Aqu ila region (~ 10- 
30 K, partially due to W40; iBontemps et alj[2010) than in Taums 
(~10-I5K). 

The mass of the Taums molecular cloud (2 — 4 x 10* M© - 
[Goldsmith et al. 2008 ; Kenvon et al. 2008) and t he Aquila molec- 
ular cloud (4 X lO* M0 - iBontemps et alj|201(lf) are comparable. 
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although the total area mapped by Herschel in Taurus is ~40% 
larger. Note that this study only presents a small part of that area. 
Aquila is also a region with on-going low-mass star formation with 
a core-mass function that peaks at ~0.7, M©. See Section l4~2l for a 
discussion of the masses of cores in this part of Taurus. 



3 COLUMN DENSITY MAP 

3.1 Pixel-by-pixel SED fitting 

The column density, A'^(H2), map of the Taurus S1-S2 region was 
reconstructed by fitting the spectral energy distribution (SED) of 
each pixel in the Herschel map. Before fitting the pixel SEDs we 
must calibrate the data in surface brightness units, account for the 
map zero-point, and then colour-correct each pixel. The default cal- 
ibration of the PACS data is already in surface brightness units. We 
converted the SPIRE data from point-source to extended calibra- 
tion an d then to surf ace brightness units using the recommended 
factors JBenddbOllh . 

Fitting pixel-by-pixel SEDs can be sensitive to systematic off- 
sets between the different wavelengths. For the SPIRE or PACS 
data alone, such offsets would not be a problem, as a common 
baseline removal/mapping technique should remove a common 
background. However, different techniques are used for the SPIRE 
and PACS instruments. Therefore, it is necessary to perform some 
measure of large-scale background correction at each wavelength, 
e.g., by masking a source and estimating a local background (e.g., 
IPeretto et al.ll2010l) . globally filtering the map, or using additional 
data to fit a zero-point. We chose the latter as it allows us to fix 
the background ove r the entire area. We therefore followed the 
method described bv lBemard et al.l (120101) to estimate the absolute 
flux level towards the S 1-S2 region from Planck and IRAS data. The 
median absolute fluxes for each wavelength are listed on line 2 of 
Tabled 

Once the backgrounds had been established, the data were 
convolved to a co mmon resolution (the SPIRE 500 /im PSF, 36" 
FWHM) using the lAniano et alj ( 1201 11) convolution kernels and 
then co-aligned on the 500 ^m pixel-grid (14" pix el width). We 
colour-coiTect each pixel using the standard SPIRE ( IBendoluOllh 
and PACS (Mulleretal. 2011) factors in an iterative cycle. The 
spectral index at each wavelength is calculated and the required 
colour correction factor is interpolated from the tabulated values. 
The correction is then applied to the data and a new set of spectral 
indices are calculated. This corrected set of indices are in turn used 
to calculate an updated correction. The spectral-index/correction- 
factor cycle is repeated until the correction factors converge or a 
maximum of 20 cycles is reached. 

For each position where there are greater than 3 a detections 
between I60-500/xm we fit the function 



iV(H2) = 



M 



mY{iiA 



(2) 



F. 



D2 



(1) 



where F,j is the flux density at frequency i/, M is the mass-per- 
pixel, Bu{T) is the Planck Function for a blackbody with temper- 
ature T, and D is the distance to the source (Hildebrand 1983). 
The dust mass opacity k^, is paramaterized as k^, oc i/^ and is ref- 
erenc ed against a value of O.I cm^ g~^ at 1 THz JBeckwith et al.l 
1 1 9901) . which assumes a gas-to-dust mass ratio of 100. The dis- 
tance to Taurus is D = 140 pc and we assume that /? = 2. There 
are therefore only two free parameters T and M for each fit. The 
mass-per-pixel is converted to a column density by. 



where A'^(H2) is the column number density of H2, my^jj, is the 
mean particle mass (mjj is the mass of a hydrogen atom and /i was 
taken to be 2.86, assuming the gas is ~ 70% H2 by mass), and A 
is the area of each pixel. 

The L evenberg-Marquar dt least-squares minimisation pack- 
age MPFIT JMarkwardll2()09b was used for the SED fitting. We use 
a single temperature fit over the range 160-500/im. This neglects 
line-of-sight temperature variations and gives a temperature that is 
representative of the entire column of material, not just its coldest 
component, seei Shettv et al.l J200 y) for a discussion. We explored a 
2-component temperature fit by including MIPS (Rieke et al. 2004|) 
lAjim data from the Taums Spitzer Legacy Program (JGiidel et al.l 
1 20071 made available from the Spitzer Science Centre). For posi- 
tions where there are 3 a detections at 24 /im, 70 /im, and 160 /im, 
it was possible to fit the SED at 24-500 /im with the sum of two 
components. In practise this only applied close to the 70 /im com- 
pact sources as these fields contain no extended-area of strong short 
wavelength emission. The cold component fit was virtually identi- 
cal to the original 1 60-500 /im fit so we only use a single compo- 
nent fit. 

Fig [2] shows the A'^(H2) and T results (for the cold compo- 
nent). Column density contours, which start at 5-(j and increase ge- 
ometrically with each increment being 1.5 times the previous con- 
tour level, are shown on the middle and right maps for reference. 
The temperature map shows minima of ~ 10 K in B18-E 2, B18- 
E I, and LI 536- 1, with the rest of the higher-density filamentary 
material lying at a temperature of '^ 12 — 14 K. Away from the fil- 
aments the temperature rises to ~ 16 K. The first analysis of Planck 
data tow ards Taurus found a similar tem p erature towards dense 
regions (Planc k Collaboration et al.l uOIlh . lOel Burgo & Laureiia 
(2005 ) used IRAS 60 pm and ISOPHOT 200 pm maps of the TMC- 
2 clump to infer two distinct dust temperature components, a warm 
20 K component which followed the 60 /^m emission and a colder 
12.5 K dense component that followed the 200 /^m emission. Their 
colder component is similar to that which we measure towards the 
filaments. 



3.2 Estimation of Uncertainty 

The intensity uncertainty of each pixel at each wavelength is taken 
as the quadratic sum of three components. The first component is 
(Jrma, the pixcl-to-pixel rms measured on the resampled, regridded 
intensity map. This was calculated by subtracting a map smoothed 
by twice the beam FWHM from the original map and measuring the 
standard deviation of the remaining intensity. The calculated values 
of (Jrms are listed in Table [T] A map uncertainty product is avail- 
able from the Herschel map making processes, but it was found 
that the PACS uncertainty maps contained striping associated with 
masked/bad pixels. While this effect is small (--^10% of the overall 
uncertainty) it is correlated with the scan direction and, if included 
in the fitting, propogates to the output column density map. There- 
fore, we instead use the per pixel rms across all wavelengths. This 
method avoids any scaling that would have been necessary due to 
confusion effects from the resampling/regridding process. 

The other two uncertainty components are related to the inten- 
sity calibration. As described above, the total intensity is the sum 
of two components ~ the Herschel measured intensity Fh at each 
pixel and the absolute background level of the map Fb ■ Each of 
these has a separate systematic fractional uncertainty. We assume 
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that the emission is extended, so the systematic fractional error Ch 
is 12% for the SPIRE bands (Valtchanov 2011) and 10% for the 
PACS 160/im aPaladini et al.ll2012l). The fractional er ror on the ab- 
solute background is Cb = 5% JBemard et al.l2010l) . 

Therefore, a pixel with total intensity F^ = Fh + Fb, will 
have a total uncertainty given by 



+ {FHCuf + {FbCbY 



(3) 



In practise the second term dominates in regions of high inten- 
sity while the last term dominates in regions of low intensity. Only 
the first term is used when determining whether the intensity of a 
pixel is statistically detected or not. The full value of a,j is used 
for the SED fitting. The uncertainty on the A'^(H2) fit is dominated 
by a systematic component of ^--^10% above a column density of 
10^^ cm~^, a value similar to the systematic error in the intensity. 
The per pixel rms of the column density, calculated in the same 
manner as for the original intensity maps, is 0.16 x 10^^ cm~^. 
The uncertainty on the T fit was found to be ~0.3K in regions 
where T < 13 K and ~1 K in regions were T > 17 K. 

We tested our fitting procedure using synthetic SEDs with 
data-points matching the Herschel bands. These were generated 
for a random-temperature Tj^ipm in the range log(5/K)-log(100/K). 
Noise was added to the data-points and the synthetic SED plus 
noise was used to obtain Tgjjgj in the above manner. In each case, 
the mean difference between the input and fitted temperature in- 
creased with the input temperature as the peak of the SED moved 
to shorter wavelengths outside of the Herschel bands. Indeed, the 
mean difference increased until a breakdown temperature, above 
which the synthetic SED could not be recovered. Our typical cal- 
ibration error is 10-12 per cent plus a component from the back- 
ground. For the conservative case of 15 per cent noise, we cannot 
recover temperatures above ~ 25 K when using just 250-500 /^m 
data-points and ~ 40 K when 160/xm data-points are added. 



4 RESULTS 

4.1 Source identification 

The maps contain many sources that are extended relative to the 
beam-size at each wavelength. The extraction of extended sources 
from regions of complex emission is a difficult process, and is one 
that many Herschel surveys will face. There has been much de- 
bate as to which method is best fo r identifying extended sources 
in sub millimetre surveys (e.g., see IPineda et al.ll2009t iReid et al.l 
l20IOf) . A comprehensive study of the current state of the art in this 
area has been carried out by the SPIRE Consortium within the Her- 
schel team, and it is still ongoing. For this first-look paper, we use 
a conservative varian t of the well-known routine CLUMPFIND 
JWilliams e t al. 1994^ which we have called the Cardiff Source- 
finding AlgoRithnQ (CS AR). CS AR also retains information per- 
taining to the structure tree of sourc es within the field mak ing 
it similar to other re cent dendrogram (IRosolowskv et al.lEOOa) or 
sub-structure codes dPeretto & Fulledl2010l) . Where CS AR differs 
is that the data are analysed as a pixel stream. This eliminates the 
need to pre-filter the map maxima and reduces the number of free 

parameters involved in its operation. 

CSAR is a seeded region-growth (c.f. lAdams & Bischoal994l) 



^ CSAR is written in IDL and can be downloaded from 

|http : \star . jmkprime ■ org] 




Fi gure 3. An illustration of how CSAR differs from CLUMPFIND (after fig 2. 
of IWilliams et aljl994h . There are four peaks A-D in this example. The red 
regions show the areas around each peak which would be assigned to them 
by the non-hierarchical version of CSAR. The four different blue shaded 
regions show the additional area around each peak which would be assigned 
to them by CLUMPFIND. The black contours show the clump divisions from 
the hierarchical version of CSAR. 



style segmentation algorithm which can be run in either hierarchi- 
cal or non-hierarchical modes. It is designed to identify contigu- 
ous regions, or clumps, on an intensity map that are considered 
significant by virtue of being: (i) resolved, namely that the clump 
area is greater than the telescope beam footprint; (ii) having no re- 
solved internal sub-structure; and (iii) having a peak-to-clump min- 
imum contrast ratio greater than 3 cr, where a is the mean pixel-to- 
pixel error. The clumps that the non-hierarchical version of CSAR 
finds are broadly equivalent to CLtJMPFIND's seed clumps before 
the ap plication of its 'friends-of-friends' algorithm (IWilliams et al.l 
[1994). The differences between CSAR and CLUMPFIND are illus- 
trated in Fig |3] which is based on fig 2. of Williams et al. (I994|) . 
The red regions are the clumps found by the non-hierarchical ver- 
sion of CSAR while the blue regions show the clumps that would 
have been found by CLUMPFIND. 

The non-hierarchical version of CSAR works by sorting each 
pixel i in the map into a list Si, in order of decreasing intensity. 
The first pixel in the list. Si , has the highest intensity on the map 
and is assigned a clump index A'^ = 1. Each pixel in the list is then 
considered in turn. The zth pixel is assigned a new clump index if 
it is a local maximum, i.e., its neighbouring pixels have not already 
been assigned a clump index. A pixel that adjoins an existing clump 
takes that clump's index, i.e., its neighbouring pixels share a single 
unique clump index. 

Clumps grow as the routine proceeds down the intensity gra- 
dient, as each new pixel is considered. This process continues until 
two clumps collide, i.e., where the ith pixel has immediate neigh- 
bours with more than one unique clump index. This point is equiv- 
alent to the position in Fig[3]where the red regions around peaks A 
and B touch. The ith pixel is then a bridging pixel between two 
clumps and the code must make a choice as to which clump it 
should be assigned to. Each clump is tested to see if it is signifi- 
cant against the three criteria described above. If both clumps are 
significant, then they are marked as finished and are considered to 
have reached their maximum size before becoming blended. They 
are then excluded from further growth. If one or both of the clumps 
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are not significant, tlien tlie clumps are merged. Every pixel is pro- 
cessed in turn until a preset minimum intensity level, or the min- 
imum on the map is reached. Each surviving unmerged clump is 
then checked to make sure it matches the resolution and contrast 
criteria and any that straddle the edges of the mapped region are 
excluded. 

The hierarchical version of the CSAR algorithm uses the un- 
blended clumps found by the n on-hierarchical version as the tips of 
a structure tree or dendrogram. IHoulahan & Scalol(ll992k proposed 
the use of structure trees to measure hierarchical structure within 
ISM clouds. These were constructed by decomposing a flux map 
into discrete grey-level masks. Clumps were found independently 
at each grey-lev el and then cross-matched with overlapping clumps 
at the N-1 level. iRosolowskv et alj J2008l) extended structure trees 
to work with three-dimensional i sosurfaces in molecular line maps. 
A true dendrogram, as shown bv lRosolowskv et al.ldlOOSi) . uses an 
arbitrary x-axis (e.g., a sorted, deprojected clump index) to show 
the hierarchical struct ure of another quantity (temperature, inten- 
sity, etc) on the y-axis JGoodman et alj2009r) . 

It is possible to use the hierarchical version of CSAR to es- 
timate the shared contribution from the surrounding emission, 
thereby overcoming any problems of the non-hierarchical version, 
which makes no attempt to deblend sources or to trace the wings 
of emission down into the surrounding cirrus. The assignment of 
clump indices and merging of non-significant clumps proceeds as 
for the non-hierarchical version except that colliding significant 
clumps are now allowed to merge. The state of those clumps is 
recorded at the point of collision. The black contours in Fig[3]show 
the clump boundaries at the point of merger. Moving down the in- 
tensity scale, clumps A and B merge at the topmost black contour, 
the combined clump containing them then merges with clump C at 
the level of second highest black contour. In this manner a network 
of merging clumps is traced. Each clump merger is a node in the 
structure tree with the unblended non-hierarchical clumps forming 
the upper-most branches of the tree. 

CSAR was run on the column density map of the SI-S2 re- 
gion shown in Figure |2] The results are shown in Figure |4] and |5] 
Figure |4] shows the dendrogram for those sources while Figure |5] 
shows a map of the sources extracted. Over-plotted on Figure[5]are 
column density contours (with levels repeated from Figure |2]l- 

CSAR found a single clump hierarchy with two dominant 
branches (northern and southern filaments), the connection be- 
tween these is shown by the lowest horizontal-line at ~ 6 x 
10^" cm^^ on Figure|4] The clumps with indices below 77 are part 
of the B18 branch while clumps with indices above 77 are part 
of the LI536 branch. For each clump, we compute a position (R.A, 
Dec) from the mean centre of the clump's half-power contour above 
its merger column density. The index of each clump is plotted at its 
position on Figure|5] The border of each clump is effectively an iso- 
contour and this level is used as that core's background level. The 
level is calculated from the ring of pixels immediately adjacent to 
each clump (the need to terminate a clump's growth with an unas- 
signed pixel automatically means that no clumps abut pixels from 
another clump). The integrated count (either mass or flux density 
depending on the map type) of each clump is then computed by 
subtracting the isocontour boundary level from each pixel in the 
clump and summing the remainder. 

The size for all CSAR sources has been taken as the decon- 
volved effective radius of a circle with the same area as the source. 
Fig|6]shows column density maps extracted from Fig|2]for a sam- 
ple four CSAR cores with approximately the same effective radius. 
The two cores on the left (cores 47 and 20) are coincident with 



the scatter of blue-triangles in Fig|7] They have masses in excess 
of 2 Mq and show evidence for being cooler at their centres (see 
Fig[2]l. The two cores on the right (cores 38 and 105) have approx- 
imately the same deconvolved size, r^ 0.06 pc, as the first two, but 
have masses less than 2 Mq . 

The two higher-mass cores have SIMBAD entries and are 
coincident with ob jects 35 and 31 from the H^'^CO"'" survey of 
lOnishi et al.l 12003). Core 47 has a Herschel mass of 5.4 Mq and 
radius of 0.06 pc while the H^"^CO^ core is it coincident with has 
a virial mass of 1.5 Mq an d a radius of 0. 04 pc. Virial masse s scale 
with R^ (c.f. lLarsonll98ll) , so scaling the'O nishi eFai] ( |2002|) virial 
mass for Core 47 to the Herschel radius gives a mass of 3.4 Mq. 
Core 20 and the H^^CO^ core it is coincident with have the same 
radius, 0.06 pc, but core 20 has a mass of 2.1 Mq and the H^^CO+ 
core has a virial mass of 0.6 A/q. Thus, the Herschel mass for 
each core is higher than the virial mass computed from H^'^CO"'' 
as would be expected if both cores were gravitationally bound. 

Conversely, the lowe r-mass core 38 is co incident with object 
31 from the CO survey of lOnishi et al.l ( 1199 6) and core 105 is not 
coincident with any extended structure in the SIMBAD databaso 
The mass and size of CO object 31 is significantly higher than 
for the Herschel core (12 versus I Mq and 0.21 versus 0.058 pc). 
Therefore, the CO core appears to be equivalent to an entire sub- 
branch of the Herschel dendrogram and not a single leaf-node. Fig- 
ure |4] and |6] show that the same Herschel observations are able to 
map simultaneously both well-known dense cores and previously 
uncatalogued lower density cores, and that the structure in the Her- 
schel maps can be characterised by CSAR. 

CSAR was applied to the column density maps in its hier- 
archical mode and structures were found on scales in the range 
0.024-2.7 pc. This result is to be expected, since molecular clouds 
generally contain structure on a variety of scales. A total of 236 
nested-sources (which we call nodes) were found in the hierarchy, 
1 15 of which do not contain resolved substructure. A dendrogram 
of the hierarchy is shown in Figure |4] 

We also ran the G ETSOURCES algorithm on our data 
JMen'shchikov et al.ll20I2r) . This routine found the same sources 
as CSAR, and others besides. GETSOURCES uses all five wave- 
lengths simultaneously to identify sources, and we saw that most of 
the additional sources it found were seen only in a subset of wave- 
lengths. Furthermore, GETSOURCES did not find the large clumps 
found by CSAR because it does not define sources as such. Any 
source containing substructure is broken into its components by 
GETSOURCES. Consequently, all of the sources discussed here were 
found by both algorithms. We believe this is the most robust way 
to proceed. 



4.2 Mass-Size relation 

Fig |7] shows the mass-size relation for the sources that we have 
found (following fig. 4 of Motte et al. 2001) . The solid red circles 
indicate the individual resolved objects with no visible substruc- 
ture and the open red circles denote the objects that we found with 
substructure. The two types of object have median masses of 0.20, 
and 25 M© respectively and median sizes of 0.055 and 0.35 pc re- 
spectively. The root clump ~ the total of everything in the hierar- 
chy - has a mass of 880 M© and an equivalent size of 2.7 pc. For 
comparison, the open and solid grey squares show data for pre stel- 
lar cores in Ophiuchus (,Motte et al.,il998.) and Orion GMotte et al.l 



2 Core 105 is coincident with the stai- 2MASS J04325387+2248375. 
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Figure 4. The full dendrogram for the S1-S2 N(H2) column density map from Fig|2] The vertical axis shows column density. Individual sources are shown by 
the vertical lines, the height of each line shows the range of column densities within that source. Horizontal lines connecting two sources together denote that 
those sources have merged at that column density. The horizontal axis is an arbitrary clump index which allows the structure tree be plotted as a dendrogram 
(i.e., a tree diagram without crossed branches). 
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Figure 5. A re-projection of the clumps as shown in the Figure |4] but now plotted over the S1-S2 N(H2) column density contours from Fig|2] The numbers 
show the arbitrary index from Figure It can be seen that, broadly, clumps with indices 1-76 are coincident with the northern filament and cores with indices 
78-H are coincident with the southern filament. The dashed contour encloses the features included in the dendrogram. 
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Figure 6. Column density maps for four CSAR cores. Thegreyscale and contours show logjQ[A'^(H2)/cm~^]. The greyscale range is 21.1-22.4. The contours 
are the same as for Fig|2] Each panel shows a source described in the text (centred at 0,0) and its sibling (the source it merges with in the dendrogram). Each 
core is labelled in red with the same index as shown in Fig|4]and its FWHM contour is shown by the red ellipse. The green-contour is the column density at 
which these sources merge. This is the same as the horizontal connecting lines on Fig|4] The blue-contour is the column density at which these merged sources 
would in turn merge with a third source. The white circle is the telescope beam FWHM. 
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Figure 7. A mass-size plot for resolved sources found in our survey, compa red to those found in other surveys using different tracers. The open and closed re d 
circles are the ne w da ta from this paper Comparison data are taken fromi lOnishi et al.l U996h ; lMotte et alj ri998l) : lMotte et alj f200ll) : IOnishi et alj J2002l) ; 
iKirk et alj )2005r) : and lSchmalzl e t al. (2010). The key at upper-left show s the meaning of each of the symbols. The sol id gray band and the dott ed line show 
the M pn oc R^^ trend for CO cores (Elmegreen & Falgaronelll996b . The soHd and fashed blue lines are fits to the lOnishi et al.l J2002h and lOnishi et alj 
119961) data respectively. The vertical red dot-dashed line is the resolution limit of the Herschel data, the diagonal red dot-dashed line is the mass sensitivity 
limit for a 3-a- Gaussian source with FWHM equal to the deconvolved radius. See text for discussion. It is assumed that the effective radius and FWHM are 
approximately equal measures of a core's size. 



UOOlh respectively, obsei-ved in the millimetre/submillimetre con- 
tinuum. Note that the same assumptions regarding mass calcula- 
tions from flux densities were used in calculating the masses for 
Orion and Ophiuch us as for the curre n t survey of Taurus . Hence, 
any offsets between lMotte et alj h998h . lMotte et alj J200lh and the 
current data are not due to different assumptions about mass opac- 
ities, or equivalent systematic uncertainties. 

We use the column density map to calculate masses. The al- 
ternative method would be to calculate the core's flux densities at 
each wavelength using a local background and then fit the SED. 



We checked what difference this second method would have on the 
leaf-node sources by fitted SEDs (using the same method as for the 
pixel-by-pixel SEDs) to background subtracted flux densities mea- 
sured using the CSAR extraction contours. It was found that using 
a local background increased the core masses by ~20-30% on av- 
erage. This increase is due to the removal of flux associated with 
warmer, overlaying material, However, using this method discards 
the information content of the dendrogram, which we wish to pre- 
serve, and requires one to have a priori knowledge of the core lo- 
cations. The advantage of the dendrogram/column density method 
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is that it preserves information on the global properties of the re- 
gion and supplies local core properties in a single self-consistent 
operation. 

Analysts of molecular cloud structure have pursued the idea 
that there exists a power-law relationship between the mass of a 
core, M, and its radius R, such that M oc R''. iLarsod ( Il98lh 
found n(H2) oc R^^'^, which is equival ent to M oc R^'^, based 
on a sa mple of mainly ^'^CO observations. lElmegreen & Falgarond 
( 119960 found a value of fc = 2.35 for a sample of clouds, 
which they studied in CO. A study of CO clouds in the inner- 

[ alaxy Molecular Ring fo und a power-law with k — 2.36 ± 0.04 
Roman-Duval et alJuOlQ ). However, the power-law exponent as 
measured by denser gas tracers (e.g. the submillimetre conti nuum) 
appears to be different. For example, Curt is & Richer! (120 101) found 
for dense cores in Perseus that a power-la w exponent v alue closer to 
~1.5 was a better fit to their data. P ineda et al.l ( 12009 ) showed that 
the exponent can be influenced by the choice of source extraction 

algorithm. 

The dataset of Elmegreen & Falgarond ( 1 199 6) and Motte et al.: 

([1998, 2001) did not include data from Taurus, so we also include 
in Fig |7] literature data from two other surveys of Taurus. The 
open b lue triangles show cores in Taurus observed by lOnishi et al.l 
( 119961) in C^*0 (1-0). These points can be fit by a power-law re- 
lationship with an exponent of k = 2.4±0.1, shown as a blue 
dashed line on Fig |7] This value of 2.4 is similar those found 
for CO clouds by other authors ( lElmegreen & Falgarond Il996t 
iRoma n-Duval et alj2O10l) . 

Onishi et alj 12003) caiTied out an H^'^CO"'" emission survey 
of Taurus. The higher critical density for collisional excitation of 
H^^CO +(1-0), of-^ 10f_cm~^ over C^*0 (1-0), of ~ 10* cm"^ 
allowed lOnishi et alj 12003) to probe denser parts of the same 
clouds. The filled blue triangles in Fig [7] show the virial masses 
calculated from these data. They also appear to follow a power-law 
trend. We fit to this trend a power-law relationship with an expo- 
nent value of k — 1.2 ±0.1 (solid blue line). Thi s value is consis- 
tent w ith that found for dense cores in Perseus bv lCurtis & Richeg 
ilOK J). The H^'^CO^ masses are virial masses (M oc R, assum- 
ing w eak or no linewidth-size correlation, c.f. fig 7. of lOnishi et al.l 
l2002h . so the value of this exponent in not unexpected. 

The lowest mass CS AR cores are consistent with the trend for 
CO clumps while the highest mass CSAR cores are in the param- 
eter regime of the H^'^CO'*' (1-0) sources. This shows that, for a 
given radius, we are sensitive to sources over a range of 2 dex in 
mass. A number of the substructured CSAR clumps are consistent 
with the k = 2.4 power-law fitted to the C^**0 sources. This means 
that we can detect the same material as the C^*0 (1-0) survey, al- 
beit over a wider range of spatial scale. Thus, the Herschel CSAR 
sources are coincident with all of the mass-size regimes shown by 
the other tracers. 



5 DISCUSSION AND CONCLUSIONS 



Ay becomes non-linear above Ay > 10 mag (Pineda et al. 201 
and b y the differing chemical history of each core (Hacar & TafalL 
UOllh . The high-spatial range of Herschel and the reliability of dust 
mass measurements makes it the perfect telescope to observe ob 
jects on a wide range of scales. 



IPineda et alJ(l2008[) show e d that C^*0 remained unsaturated below 
^„~10 and lPineda et al.l 12010! ) showed that CO and dust masses 
for parsec sized clumps were in good agreement. We see a similar 
effect in Fig |7] where the clumps from our dendrogram that have 
radii of ~0.3 pc have masses which are comparable to CO clumps 
of the same size lOnishi et al.l 19961) . However, the case is different 
at smaller scales where interpretation of individual cores is ham- 
pered by the spatial filte ring of ground-based bolometer observa- 
tions IPineda et alj|201(lh , the relationship between CO mass and 



iKauffmann et alj 120 ICl ) produced a mass-size diagram using 
a dendrogram technique for multiple clouds. They concluded that 
the mass-size data were not consistent with a constant column den- 
sity relationship (k = 2) and that cores in low-mass star formation 
regions had masses that lie below a threshhold with an exponent 
of fc = 1.33 (close to the solid blue line in Figure |7j. However, 
IKauffmann et alj 12010 ) had to use two different data-sets to bridge 
the same parameter space. What Herschel adds to this type of anal- 
ysis is its ability to sample a larger range of spatial scales in a sin- 
gle map and to extract a greater dynamic mass range of object for a 
given radius. Nevertheless, the linear part of the mass-size relation 
in the top-left of Figure [T] as the structure tree folds into the root 
node, does appear to be inconsisten t with a global value of fc = 2 
in agreement with lKauffmann et alj 120101) . 



As mentioned previouslv. iMotte et alj ( 1200 II) found that cores 
seen in the submillimetre continuum (squares on Fig|7) and clumps 
seen by molecular line tracers such as CO (grey shaded area on 
Fig IT} lie in very different areas of the mass-size plane. Our Her- 
schel data allows is to study the two possible explanations for this 
difference: either selection effects are at play, or there are two fun- 
damentally different populations of object , i.e., low density non- 
star forming starless cores (c.f . Ward-Thompson et al.lEOlOh. an d 
higher density prestellar cores dWard-Thompson et al.ll994l . l2007h . 



If the latter explanation were correct, then the populations 
would remain distinct in the mass-size plane, even after selection 
effects were corrected. However, Herschel is sensitive to much 
lower column density material than previous ground-based sub- 
millimetre continuum arrays. Consequently, the sensitivity limit of 
Herschel continuum data lies further to the lower right on Fig |7] 
Hence, Herschel effectively removes the old selection effect prob- 
lems, and so its data should be able to settle this issue. 

Previous results from Herschel observations of other regions 
have been moderately inconclusive in resol ving this issue. In more 
activ e star-forming regions such as Aquila ( lAndre et al.ll2010) and 
Vela IGiannini et alj|2012r) . the cores found by Herschel lie in the 
same part of the diagram as the ground-based submillimetre con- 
tinuum cores, indicating that they are mostly prestellar cores. In 
non star-forming regions, such as the Polaris, the cores found by 
Herschel lie in the same band of the diagram as the CO clumps 
dAndre et alj|2010|). indicating tha t they were mostly CO clumps 
(c.f. Ward-Thompson et al. 2010). A gap remained in the mass- 
size relation plot between these two samples. 

Now we can see from Fig|7]that the Herschel cores in the Tau- 
rus star-forming region lie largely in the gap between the ground- 
based continuum cores and CO clumps. Their appearance rules out 
the possibility that there are two distinct populations of cores. In- 
stead, the Herschel data indicate that all cores are part of the same 
population. Furthermore, this conclusion leads to the hypothesis 
that lower density cores can evolve into prestellar cores by mi- 
grating across this diagram before collapsing to form protostars. 
A suggested form of th is evolution has been previously proposed 
by Si mpson et alj ( 120111) . Further results from the Herschel Gould 
Belt Survey in other regions should help to resolve the exact man- 
ner of the evolution of prestellar cores. 
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